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Abstract 

. . . 

Given two triangulations of a convex polygon, computing the minimum number of 

flips required to transform one to the other is a long-standing open problem. It is not 
known whether the problem is in P or NP-complete. We prove that two natural gen- 
eralizations of the problem are NP-complete, namely computing the minimum number 
of flips between two triangulations of (1) a polygon with holes; (2) a set of points in 
^ the plane. 

u 

g 1 Introduction 

Given a triangulation in the plane, a flip operates on two triangles that share an edge 
and form a convex quadrilateral. The flip replaces the diagonal of the convex quadrilateral 
by the other diagonal to form two new triangles. A sequence of flips can transform any 
triangulation to any other triangulation — this is true for triangulations of a convex polygon, 
and more generally for triangulations of a polygonal region with holes, which includes the 

^ case of triangulations of a point set. 

In this paper we investigate the complexity of computing the flip distance which is the 
minimum number of flips to transform one triangulation to another. This is particularly 
interesting for convex polygons, where the flip distance is the rotation distance between two 

^ binary trees (see below). 

The main result of our paper is that it is NP-complete to compute the flip distance 
between two triangulations of a a polygon with holes, or a set of points in the plane. 



1.1 Flip distance and rotation distance 

Balanced binary search trees are a widely used data structure. One way to make a rooted 
binary search tree balanced is using an operation called a rotation [5J. Despite being very 
simple and fundamental, the rotation operation is not completely understood. In particular, 
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the complexity of the problem of computing the minimum number of rotations needed to 
convert one rooted and labelled binary search tree to another, called the rotation-distance, 
has been open since 1987 [HIE]. It is not known if the problem is NP-complete. The problem 
is also closely related to Sleator and Tarjan's famous dynamic optimality conjecture [18]. 

Rooted binary trees are of interest to researchers in the field of bioinformatics as well. An 
evolutionary tree is a rooted binary tree and constructing an evolutionary tree that fits a given 
set of species given some information about their DNA sequences is a widely studied problem. 
Evaluating the effectiveness of methods of constructing evolutionary trees leads naturally to 
the problem of measuring how similar two trees are to each other. Several measures of 
similarity have been used including the rotation distance between the two trees [6]. 

There is a bijection between binary trees with n — 1 labeled leaves and triangulations of 
an n-vertex convex polygon. Moreover, a rotation in the tree corresponds to a flip in the 
polygon. Thus, computing the rotation distance between two trees is exactly equivalent to 
computing the flip distance between two triangulations of a convex polygon. See [19]. 

1.2 Generalizations and related work 

Flip distance between triangulations of a convex polygon and rotation distance between 
binary trees have been well studied in the past. Several results deal with the combinatorics 
of the flip operation. Sleator et al. [19] proved that for large values of n, the flip distance 
between two triangulations of an n-gon is at most 2n — 10 and that occasionally 2n— 10 flips 
are necessary. Other results deal with the flip graph, i.e., the graph where nodes correspond 
to triangulations of an n-gon and an edge between two nodes denotes the fact that the 
corresponding triangulations are one flip apart. For example, Lucas p2] showed that the flip 
graph is hamiltonian. See also Eppstein [8]. 

Flips have been studied in more general settings as well. Dyn et al. [7] proved that any 
two triangulations of a simple polygon can be transformed into one another using flips. They 
proved that the same holds even for two triangulations of a simple polygon with points inside 
it. Lawson [TO] proved an upper bound of 0{n 2 ) flips needed in any such flip sequence. 
Hurtado et al. [9] proved that the bound is tight asymptotically. 

Triangulations of polygons with interior points have been further generalized to triangu- 
lations of simple polygons with polygonal holes, also called polygonal regions. It is known 
that two triangulations of the same polygonal region can be transformed into each other us- 
ing flips (see [13]). Note that a one- vertex polygonal hole is just a point. Thus triangulations 
of polygons with interior points are a special case. 

Flips have also been studied in a more combinatorial setting. For example, given a 
maximal planar graph, we can define a flip as replacing an edge with another so that the 
resulting graph is also maximal planar. Wagner [20] proved that given two maximal planar 
graphs G\ and G2, there always exists a sequence of edge flips that transforms G\ into a 
graph isomorphic to G2- Combinatorial bounds on the number of flips required have also 
been studied and the best known upper bound is by Bose et al. [3] of 5.2n — 24.4. 

In the combinatorial setting, we have the choice of labelled vs. unlabelled graphs. Sleator 
et al. [TB] proved that O(nlogn) flips are sufficient to transform one labelled maximal planar 
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graph with n vertices into another with the same vertices, and Q(n\ogn) flips are sometimes 
necessary. Needless to say, this is a huge area with numerous directions of investigation. 
Bose and Hurtado [2] provide a survey. 

Regarding the question of actually computing the flip distance, to the best of our knowl- 
edge, only triangulations of convex polygons have been studied and the question has been 
open since 1987 [T71 The best known result is a trivial factor-2 approximation algorithm, 
which can be improved under certain assumptions regarding the input [TT] . Recently it was 
proved that the problem is fixed parameter tractable in the flip distance [I]. No hardness 
results are known either. 



2 Triangulations of polygonal regions 

Theorem 1. The following problem is NP-complete: Given two triangulations of a polygonal 
region with holes and a number k, is the flip distance between the two triangulations at most 
k? 



2.1 Proof idea 

Note that the problem lies in NP. We prove hardness by giving a polynomial time reduction 
from vertex cover on 3-connected cubic planar graphs, which is known to be NP-complete [H 



The idea is to take a planar straight-line drawing of the graph and create a polygonal 
region by replacing each edge by a "channel" and each vertex by a "vertex gadget" . We then 
construct two triangulations of the polygonal region that differ on the channels, and show 
that a short flip sequence corresponds to a small vertex cover in the original graph. 

We begin by describing channels and their triangulations, because this gives the intuition 
for the proof. A channel is a polygon that consists of two 7-vertex reflex chains joined by 



two end edges, as shown in Figures 1(a) and |l(b) . Note that every vertex on the upper 



reflex chain sees every vertex on the lower reflex chain and vice versa. We identify two 



triangulations of a channel: a left-inclined triangulation as shown in Figure 1(a); and a 



right-inclined triangulation as shown in Figure 1(b) 

A channel is the special case n = 7 of the polygons H n of Hurtado et al. [9]. They prove in 
Theorem 3.8 that the flip distance between the right-inclined and left-inclined triangulations 
of H n is (n — l) 2 . We include a different proof in order to generalize: 

Property 1. Transforming a left-inclined triangulation of a channel to a right-inclined tri- 
angulation takes at least 36 flips. 

Proof. In any triangulation of a channel, each edge of the upper reflex chain is in a triangle 
whose apex lies on the bottom reflex chain. This apex must move from lower right {Bj) 
to lower left (-Bi), in order to transform the left-inclined triangulation to the right-inclined 
triangulation. Similarly, each edge of the lower reflex chain is in a triangle whose apex lies 
on the upper reflex chain, and must move from upper left to upper right. However, one flip 
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can only involve one edge of the upper chain and one edge of the lower chain (no other 4 
vertices form a convex quadrilateral), and thus can only move one upper and one lower apex, 
and only by one vertex along the chain. Twelve triangles times six apex moves per triangle 
divided by two apex moves per flip gives a lower bound of 36 flips. □ 

We now show that the number of flips goes down if a channel has a cap, an extra vertex 
that is visible to all the channel vertices, as shown in Figure 1(c). 

Property 2. The flip distance from a left-inclined to a right-inclined triangulation of a 
capped channel is 24- 



Proof. The "canonical" triangulation shown in Figure [1(d) is 12 flips away from both the left- 
inclined and the right-inclined triangulations of a capped channel: To flip the left-inclined 
triangulation to the canonical triangulation, flip edges A\B\, . . . , AiB? followed by edges 
A2B7, . . . , AqBj in that order. Similarly for the right-inclined triangulation. 

For the lower bound, we follow the same idea as above. In any triangulation, each edge 
of the upper [lower] reflex chain is in a triangle whose apex is either the cap or a vertex of 
the lower [upper] chain. There are only two kinds of flips: (1) a flip involving the cap vertex, 
an edge of one chain, and a vertex of the other chain; and (2) a flip involving one edge of 
each chain. A flip of type (1) moves the apex of only one triangle, and moves the apex to 
or from the cap. If a triangle is altered by a flip of type (1) then at least two such flips are 
required, one to move the apex to the cap and one to move the apex from the cap. If a 
triangle is only altered by flips of type (2), then, as above, it costs 3 flips to get the apex to 
its destination. Thus the 12 triangles require at least 24 flips. □ 

With the description of channels in place, we now elaborate on the idea of our reduction. 
We create a polygonal region by replacing each edge in the planar drawing by a channel, 
and each vertex by a vertex gadget. We make two triangulations of the polygonal region. 
In triangulation T\ all edge channels are left-inclined and in T2 all edge channels are right- 
inclined. The triangulations are otherwise identical. We design vertex gadgets so that making 
a few flips in a vertex gadget creates a cap for a channel connected to it. Since transforming 
a channel from left-inclined to right-inclined is less costly if it is capped, the minimum flip 
sequence that transforms all the channels is obtained by choosing the smallest set of vertices 
that covers all the edges and using them to cap all the channels. Thus, intuitively, a minimum 
flip sequence corresponds to a minimum vertex cover. 

One complication is that we cannot construct a vertex gadget for a sharp vertex — a 
vertex of degree 3 where one of the three incident angles in the planar drawing is > 71. 
Therefore, we first show how to eliminate sharp vertices. Let G be our given 3-connected 
cubic planar graph. Using a result of Rote [Hj, we can find, in polynomial time, a strictly 
convex drawing of G on a polynomial-sized grid. (In fact, Tutte's algorithm would also suffice 
for our purposes.) Strictly convex means that each face is a strictly convex polygon. Thus 
the only sharp vertices of this drawing are the vertices of the outer face. We replace each 
sharp vertex v by a 3- vertex chain v 1 , v 2 , v 3 as shown in Figure [2] We claim that G has a 
vertex cover of size < k if and only if the modified graph has a vertex cover of size < k + 1, 
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(a) A left-inclined channel. (b) A right-inclined channel. 




(c) A capped channel. (d) The canonical channel 




(e) Narrow and wide mouths. 
Figure 1: Channels 
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where t is the number of vertices on the outer face of G. This is because there are only two 
ways to cover the edges (vi,v 2 ) and (^2,^3): {vi,Vs}, which corresponds to v being in the 
vertex cover of G; or {^2}, which corresponds to v not being in the vertex cover of G. 




Figure 2: Eliminating sharp vertices 



2.2 Details of the reduction 



For the remainder of the proof we will assume that we have a graph G with vertices of degree 
2 and 3, and a straight-line planar drawing F of the graph on a polynomial sized grid with 
no sharp vertices. 

For each channel, we define its narrow mouth (the shaded region) and wide mouth (the 
dotted lines) as shown in Figure 1(e). Any point lying inside the narrow mouth and outside 



the channel can be a potential cap for the channel. We argue below that a vertex outside 
the wide mouth does not help reduce the flip distance. 

Each of 



We now describe the triangulated vertex gadgets. See Figures 3(a) and 3(b) 



the 2 or 3 channels attached to the vertex gadget will have one potential cap. We place a 
convex quadrilateral CDEF with diagonal CE, called the lock, that separates each channel 
from its potential cap. Thus the lock CE must be flipped, or "unlocked", in order to cap 
any channel. 

For the degree-2 gadget (see Figure 3(a) ), place point C in the smaller angular sector (of 



angle < 7r) between the two channels, so that C is outside the wide mouths of both channels. 
Place points D, E, and F in the other angular sector, with D inside channel l's narrow 
mouth and outside channel 2's wide mouth, E outside the wide mouth of both channels, 
and F inside channel 2's narrow mouth and outside channel l's wide mouth. Triangulate as 
shown. Thus D is a potential cap for channel 1 and F is a potential cap for channel 2. 
For the degree-3 gadget (see Figure 3(b)), note that because the vertex is not sharp, 



the mouth of each channel exits between the other two channels. We place vertices in the 
angular sectors as shown in the figure. Place D inside the intersection of the narrow mouths 
of channels 1 and 2, and outside the wide mouth of channel 3. Place F inside channel 3's 
narrow mouth and outside channel 1 and 2's wide mouths. Place C and E outside the wide 
mouths of all the channels and triangulate as shown. Thus D is a potential cap for both 
channel 1 and 2 and F is a potential cap for channel 3. 

Observe that every channel is blocked from its unique potential cap by exactly 3 edges. 



(For example, in Figure 3(b) , channel 1 is separated from its potential cap D by edges FA, 
FE, and CE.) Observe furthermore that for each vertex gadget, the sets of blocking edges 
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of the channels have one edge in common, namely the locking edge CE, and are otherwise 
disjoint. These properties are crucial for correctness. 

We will say that a vertex gadget is locked if the diagonal CE exists and unlocked other- 
wise. We first show what is possible with unlocked vertex gadgets. 

Property 3. If we unlock a vertex gadget then, for each channel attached to it, there is a 
sequence of 28 flips that transforms the channel triangulation and returns the vertex gadget 
to its (unlocked) state. 

Proof. We first claim that there is a 2-flip sequence that caps the channel. We enumerate 
the possibilities (refer to Figure [3]). Note that we handle channels one at a time, not simul- 
taneously. For the degree-2 gadget: flip CF followed by CA for channel 1; flip CD followed 
by CB' for channel 2. For the degree-3 gadget: flip EE followed by FA for channel 1; flip 
CF followed by CA' for channel 2; flip ED followed by EA" for channel 3. Once the channel 
is capped, we can transform the left-inclined triangulation to the right-inclined triangulation 
in 24 flips by Property [2j Then we undo the 2 flips that capped the channel. The total 
number of flips is 28. □ 

Next we show some lower bounds on the number of flips. First we note that the proof of 
Property [T] carries over to the following: 

Property 4. Transforming a left-inclined triangulation of a channel to a right-inclined tri- 
angulation takes at least 36 flips even in the presence of other vertices, so long as the other 
vertices lie outside the wide mouths at either end of the channel. 

We now consider what happens when we unlock some vertex gadgets. Let T[ be the 
triangulation obtained from T\ by unlocking some vertex gadgets. Let T' 2 be the triangulation 
obtained from T2 by unlocking the same vertex gadgets. Let C be the set of channels that 
have a locked vertex gadget at both ends. Then: 

Property 5. The number of flips required to transform T[ to T' 2 is at least 28\E — C\ + 36|C| . 

Proof. Consider a channel of C, with a locked vertex gadget at both ends. The cap vertices 
of the channel are not useable. By construction, the other vertices are outside the wide 
mouths of the channel. Therefore, by Property |4j we need 36 flips to transform it. 

Consider the channels with an unlocked vertex gadget at one end. We only save flips by 
capping the channel. To do this, we must flip the two edges that block the channel from its 
cap. Because the edges that block one channel are disjoint from the edges that block any 
other channel, we must do two flips per channel, and we must re-flip those edges to return 
to the original state. Finally, by Property [2] it takes at least 24 flips to transform a capped 
channel. (Note that the proof of Property [2] carries over even if the channel is capped at 
both ends.) The total number of flips is 28 per channel. □ 
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Channel 2 



Channel 2 




(a) Degree 2 (b) Degree 3 



Figure 3: Gadgets for vertices 
2.3 Putting it all together 

Lemma 1. G has a vertex cover of size < k if and only if the flip distance between the two 
triangulations T\ and T 2 is < 2k + 28|i£|. 

Proof. Suppose that G has a vertex cover of size k. Unlock the corresponding k vertex 
gadgets. Each edge channel has an unlocked gadget at one end, so by Property [3] we can 
transform between the two triangulations of the channel in 28 flips. When all channels have 
been transformed, we relock the k vertex gadgets. The total number of flips is 2k + 28|i?|. 

For the other direction, suppose that there is a flip sequence between 7\ and T 2 of length 
< 2k + 28|i£|. Let L be the set of vertices whose gadgets are unlocked in the flip sequence. 
Let C be the set of edges not covered by vertex set L. By adding one vertex to cover each 
edge of C, we observe that G has a vertex cover of size \L\ + \C\. Thus it suffices to show 
that |L| + |C| < k. By Property [5] the number of flips is at least 2|L| + 36|C| +28( | E-C\) > 
2\L\ + 28\E\ + 8|C|. By assumption, the number of flips was < 2k + 28\E\. Therefore 
2\L\ + 8|C| < 2k, which implies that \L\ + \C\ < k, as required. □ 

The last ingredient of the NP-completeness proof is to show that the reduction takes 
polynomial time. We need the following claim. 

Claim 1. The size of the coordinates used in the construction is bounded by a polynomial in 
n. 
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Proof. We give the main idea here, with further details in the appendix. We begin with a 
straight line drawing of a graph on a polynomial size grid. Expand the grid, and allocate a 
square region around each vertex for the vertex gadget. Expand each edge to two parallel 
line segments. These line segments will become the channel, but for now, the reflex vertices 
of the channel are all collinear, which means that the channel's wide mouth is equal to its 
narrow mouth. The points C, D, E, F of the vertex gadget go in feasible regions defined by 
the wide and narrow mouths (e.g. in the 3-channel gadget, point D lies in the narrow mouth 
of channels 1 and 2, but outside the wide mouth of channel 3). We make the channels narrow 
enough so that all the feasible regions intersect the region allocated to the gadget. We claim 
that we can choose the channel end points A, B, A', B', A", B" on the expanded grid so that 
the resulting channels satisfy this property. 

Now we pick points C, D, E, F inside the appropriate regions. Because the boundaries 
of the feasible regions are determined by pairs of points on the expanded grid, the new 
points can be chosen to have polynomial size (because solutions to linear programs have 
polynomial size as shown in Theorem 10.1 of [H]). Finally we place the reflex points of each 
channel. The feasible region wherein each set of reflex points can be placed is bounded by 
lines through pairs of points already placed. Thus, we can choose reflex points of polynomial 
size. □ 



3 Triangulations of point-sets 

We prove the NP-hardness of computing the flip distance between two triangulations of 
a point set by reducing from computing the flip distance between two triangulations of a 
polygonal region. Given two triangulations T\ and T 2 of a polygonal region R, we triangulate 
all the holes and pockets of R the same way in both triangulations. Next, we repeat each 
edge on the boundary of the holes and pockets many times (as shown in Figure [1]) so that 
dismantling a boundary edge requires a large number of flips. This gives two triangulations 
of a point set such that the flip distance between the two triangulations is the same as the 
flip distance between the original T\ and T 2 . 

Theorem 2. The following problem is NP-complete: Given two triangulations of a point set 
in the plane, and a number k, is the flip distance between the triangulations at most k ? 



4 Conclusion 

We have shown that it is NP-complete to compute the flip distance for triangulations of a 
polygonal region, or a point set. The problem remains open for a convex polygon, or a simple 
polygon, and also for more combinatorial objects such as labelled and unlabelled maximal 
planar graphs. 

Acknowledgements. We thank Therese Biedl for helpful suggestions. 
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Figure 4: Repeating edges on the boundary of holes. 
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A Proof of Claim H] 



As mentioned in Section 2.2, our proof consists of three steps: 



1. Draw the channels with straight and parallel chains so that the feasible regions are 
non-empty and contained in a small square surrounding the original vertex. Moreover, 
the points A, B, A', B', A", B" get polynomial sized coordinates. 

2. Obtain a point with polynomial sized coordinates inside each feasible region. 

3. Make all the chains reflex. 

For the first part, we begin by placing a square of side c around each vertex v with v 
at its center (Figure [5J, such that c is at least a constant factor smaller than the smallest 
edge in the drawing T and can be written with polynomial number of bits. If an edge passes 
through a corner of the square, then we slightly increase one of its sides. Our aim is to 
find the points A, B, A', B', A", B" on the boundary of the square. (Note that even though 
the edge AB will not be orthogonal to the two chains of the channel, the properties of the 



channel that we proved in Sections 2.1 and 2.2 will still hold. 



The edges and their extensions (the dotted lines in Figure [5J) intersect the square at 
points whose coordinates are polynomial size. Let S be the set of intersection points and 
the corners of the square. For the edge corresponding to channel 1, consider the point p\ 
where it intersects the square and find the point p other than itself in 5* that lies on the same 
edge of the square and is closest to it. Setting A to be the point on the boundary of the 
square a distance ppi/3 away from pi towards p and B the symmetric point on the opposite 
side determines the channel and its width. Do the same thing at the other end of the edge 
corresponding to channel 1 and obtain another width. Finally, pick the narrower of the two 
options for channel 1. Since A and B lie on the edge of the square and their distance to p\ 
is polynomial, we need polynomial number of bits to express the coordinates of A and B as 
well. Repeat the above for A',B',A" and B" . 

Since all the possible intersection points between the upper and lower chains of the 
channels occur inside the square, all the feasible regions have non-empty intersections with 
the interior of the square. 

Now, since the feasible regions are non-empty and are defined by linear inequalities with 
polynomial sized coefficients, using the theory of linear programming, we find a point with 
polynomial size coordinates inside each feasible region. Finally, to make the chains reflex, we 
find a location for each point on it one by one. Each point on the chain has a feasible region 
now defined by two kinds of constraints: 1) for every point outside the channel, if the point 
was inside the narrow mouth, it should remain inside and if it was outside the wide mouth, 
it should remain outside and 2) the new location of the point should maintain the reflexivity 
of the chain. These two constraints are also linear and thus we can find polynomial size 
coordinates for each point on the chains. This concludes our proof. 
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Channel 2 



Channel 2 




Channel 1 




Channel 1 



Channel 3 



(a) Degree-2 (b) Degree-3 

Figure 5: Constraints for vertex gadgets. 
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